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SUMMARY 


The meteoroid environment is of considerable scienti- 
fic and aerospace interest. With the current interest in manned 
interplanetary missions with aphelion inside the asteroidal 
belt, an estimation of the deep space meteoroid environment is 
necessary. Inasmuch as inadequate information exists regarding 
the number density of meteoroids in the asteroidal belt, this 
study was undertaken in an effort to make a constribut ion in 
this direction. 

More specifically, the physical significance of popu- 
lation index type particle distributions is examined. These are 
distributions where the number density f(m)dm of particles per 
unit volume with a mass in the range m to m + dm is given by 


f(m)dm = Am a dm 


(S-l) 


where A and a are constants, the latter known as the population 
index. 


Because of its simplicity, eq. S-l is very useful 
since if f(m) is known at two mass values, then both A and a are 
defined and extrapolation is possible. In the absence, however, 
of a priori physical reason for the existence of populations of 
the form eq. S-l, such extrapolations are questionable at best. 
The purpose of the present study is to discuss a particular cir- 
cumstance under which a population of meteoroids is correctly 
described by eq. S-l, and to compare the results with available 
experiment . 

The following equation constitutes the formulation of 
the physical model considered, and will be called the collision 
equation: 



Bellcomm, Inc. 

Summary ( cont ? d ) 


2 


time rate of change in the number density of 
particles having a mass in the range m to 
m + dm 


number of particles removed (per unit time) 
from the mass range m to m + dm due to the 
erosive influence of collisions with compara- 
tively small particles* 


number of particles removed (per unit time) 
from the mass range m to m + dm due to ca- 
tastrophic collisions (these result in the 
complete disruption of the "test" objects 
with masses in the range m to m + dm) 


(S-2) 


+ 


number of particles, in the mass range m to 
m + dm, created (per unit time) due to colli- 
sional fragmentation between objects with 
sufficiently large masses 


It is assumed that the population has reached a 
steady-state value, i.e., the time derivative term on the left 
hand side of eq. S-2 is zero, or very small. It is then shown, 
in the text, that a solution of the form eq. S-l satisfies the 
collision equation provided that the population index a is 
given by 


1.75 < a < 2 


(S-3) 


in the mass range 


uT' << m << AM /I” 

oo 


(S-4) 


*This term also includes the converse effect of slightly 
larger particles eroding into the mass range m to m + dm. 






Bellcomm, Inc. 

Summary (cont’d) 


- o - 


where y is the mass of the smallest object not blown away by 
radiation pressure, yT' is the largest mass completely shatter- 
ed by a projectile with mass y and AM^/r’ is the mass of the 

largest fragment when the largest mass in the sample, M ro , is 
shattered catastrophically. 

The collision equation eq. S-2 cannot be satisfied by 
trial solutions of the form eq. S-l for values of a beyond the 
limits given by eq. S-3. For a > 2 erosion by the smallest par- 
ticles dominates and for a < 1.75 collision products from the 
most massive objects cause evolution of the mass distribution 
with time. 


Comparison of this result with observational infor- 
mation regarding the near earth environment is given in Fig. 1. 
It can be seen, from this figure, that the gross features of 
the distribution are reproduced by the present model. 

The distribution of asteroids catalogued by Kuiper 
et al is given in Fig. 2. A least squares fit to the number 
density of the observed asteroids gives 


a = 1.80 + .04 


(S-5) 


in good agreement with our result, eq. S-3- Comparison with the 
distribution of lunar craters and meteorites is also found to be 
favorable . 


It therefore appears that an extrapolation of the num- 
ber density of asteroids into mass ranges of much smaller object 
is a reasonable first approximation to their distribution. This 
extrapolation is sketched in the summary chart Fig. S-l. It can 
be seen, from the figure, that the number density of asteroidal 
debris as obtained by the present study is considerably lower 
than the NAA engineering model. 



LOG CUMULATIVE FLUX IN METERS ' 2 SEC' 





FIGURE S-l 
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ABSTRACT 


A collisional model of meteoroids is formulated. An 
equation is derived which describes the evolution of a system 
of particles under the processes of particle destruction due 
to collisions and particle creation due to fragmentation dur- 
ing collisions. If the system of particles is assumed to have 
reached steady-state conditions it is found that a particle 

number density function of the simple form Am -ot dm (m is parti- 
cle mass, A and a are constants) satisfies the equation pro- 
vided that a has a value in the range ^-(n + 5/3) < a < 2 where 

n is a material parameter having a value of about 1.8. Compar- 
ison with near earth meteoroid fluxes indicates that the gross 
features of the observed distribution are reproduced by the 
present model. Agreement is also found with the distribution 
of the catalogued asteroids. 
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COLLISIONAL MODEL OF METEOROIDS 


1.0 INTRODUCTION 


The distribution of interplanetary debris is of con- 
siderable scientific and engineering interest. As the smallest 
members of the solar system, these objects are of interest to 
the astronomer and since they can collide with a spacecraft, 
they are of equally strong interest to the aerospace engineer. 
With the current interest in a manned Mars flyby mission which 
is planned to enter the asteroidal belt at aphelion, the esti- 
mation of the distribution of debris in the asteroidal belt is 
one of the most important environmental problems to be solved. 
This paper treats theoretically the collective dynamical pro- 
perties of a population of debris particles in orbit around 
the sun undergoing mutual collisions and subsequent fragmenta- 
tion; a physical basis for estimating particle distributions is 
therefore established. Application of the present model to the 
flux of meteoroids near earth serves as an independent check; 
the gross features of the particle flux into the earth’s atmo- 
sphere are found to be reproduced with the present model. 

Section 2 of this paper is a discussion of the frag- 
mentation of rock during hypervelocity impact. A mathematical 
formulation of the fragmentation process, based on results by 
Gault, Shoemaker and Moore (1963) is developed. 

The collision equation governing the dynamics of a 
collection of particles undergoing mutual collisions and frag- 
mentation is developed in Section 3. This is an integro dif- 
ferential equation for the number density function of particles 
having a mass in the range of m to m + dm. 

In Section 4, a trial solution of a population index 

form 


f(m)dm = Am a dm 


( 1 ) 


where ffc)om is the number density function of particles per unit 
volume having a mass in the range m to m + dm and a is the popu 
lation index, is then substituted into the collision equation. 
It is found that a solution of this form (eq. 1) solves the col 
lision equation only if the latter has reached a steady-state 
condition in time, i.e., when the particle removal rate equals 
the particle creation rate. 
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In Section 5 S th.e results of Section 4 are compared 
with the distribution of near earth meteoroids of presumably 
cometary origin. It is found that the gross features of the 
observed distribution are reproduced with the present model. 

Section 6 is an application of the model for the esti- 
mation of the number density of debris in the asteroidal belt. 

In Section 7 agreement is sought and found between the solution 
of the collision equation and the distribution of known aster- 
oids. Comparison with the distribution of lunar craters and of 
meteorites is also found to be favorable. 

2.0 THE CRUSHING LAW 


Interplanetary space contains a very large number of 
objects having different masses and orbits. These objects are 
believed to frequently collide with each other inelastically . 
When such a collision occurs at a sufficiently high relative 
velocity, fragmentation results. In the present study, the 
relative velocities will be comparable to those of dust parti- 
cles in space traveling in different but intersecting orbits. 
This means that the impact velocity will be of the order of 
kilometers per second and, hence, sufficiently high to cause 
fragmentation . 

Regarding the mass distribution of fragments produced 
during impact, the following type of crushing will be assumed: 


g(m; MjIV^dm = C(M,M 2 )m n dm 


( 2 ) 


Here, g(m; is the number of particles having a mass 

between m and m + dm produced during the impact of a mass M 
with another, larger mass M 2 . The coefficient C(M,M 2 ) is a 

function of the colliding masses and n is a constant. 

We first consider the case when the target mass is 
very great compared with the projectile particle, i.e., 

(M 2 /M) -*■ «. This is the frequently considered problem of high 

velocity impact into a semi-infinite target. The quantity g is 
now a function of M and m only and eq. 2 becomes 


g(m; M,°°)dm = C(M,®)m 11 dm 


( 2 * ) 


This particular crushing law for a semi-infinite target is 
based on experiment (Gault, Shoemaker and Moore, 1963) and 
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observation over a limited number of cases. Use of a particular 
crushing law is one of the major assumptions in this paper. 
However, since evidence supports a crushing law of the general 
form of equation 2* during hypervelocity impact, it will be 
adopted here to estimate the distribution of particles resulting 
from inelastic collisions at orbital velocities. 

In order to estimate C(M,°°) the following integral has 
to be evaluated 


M 


e 



m g(m; 


M , 00 ) dm 




y 


(3) 


where M g is the total ejected mass, is the largest and y the 

smallest fragment produced; g(m; M,<=°) is given by eq. 2. Sub- 
stitution of eq. 2* into eq. 3 and subsequent integration gives 


C(M,°°) 


(2-n)M e 



, 2 -" 


(4) 


Gault, Shoemaker and Moore (1963) found that n is about 1.8; y 
is a submicron particle (Gault and Heitowit, 1963) correspond- 

_l4 

ing to a mass of the order of 10 kg or smaller. In the pre- 
sent treatment, the particle masses in the range of interest 

-14 

will be very much larger than 10 kg; one therefore obtains 
(since n < 2) 


.. 2— ri 

M, >> 
b 


2-n 


( 5 ) 


The denominator in equation 4 can, therefore, be 
simplified and one obtains 
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Both the total ejected mass M g and the largest frag- 
ment are functions of the mass M of the colliding particle 
and will be assumed to have the form 

M = TM << 
e 2 

(7) 

M, = AM << M 0 
b 2 


where r and A are both functions of the impact velocity and 
material properties of the target as well as the projectile 
but not of their masses. is the mass impacted by M and the 

double inequality sign reflects the fact that the coefficients 
r and A refer to a semi-infinite target. It can be seen from 
eq. 7 that r is the total ejected mass per unit projectile 
mass and A is the mass of the largest fragment per unit pro- 
jectile mass. 

The use of eq. 7 is based on results from hypervelo- 
city experiments discussed by Gault et al (1963). These 
authors find that the total ejected mass as well as the mass 
of the largest fragment during hypervelocity cratering into 
basalt is proportional to the projectile kinetic energy, and 
hence, to the projectile mass. These experiments were conducted 
at impact velocities over a range not exceeding 10 Km/sec and 
over a range of projectile kinetic energies from 10 joules to 
4 

10 joules, approximately. 

If 


rM £ m 2 


then eq. 7 breaks down because hypervelocity impact into a 
relatively small target differs from the former (.semi-infinite 
target) situation since the shock formed during impact will be 
reflected back toward the impact area rather than propagated 
away to infinity (i.e., dissipated). This is particularly 
significant for stones fracturing easily under tension. For 
these objects, a mass 


M 2 >> rM 
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can still be completely shattered by the shock wave (generated 
during the event) which is reflected at the free surfaces and 
propagated inward as a tension wave. 

In the absence of sufficient factual information de- 
scribing this catastrophic process, the following will be 
assumed: 

(i) the largest mass M 2 completely shattered by M 
is given by 

m 2 = r'M 

with M 5 M + M_ and r ' > r 
e 2 

(ii) when 


m 2 > r'M 


the semi-infinite target relations are valid. 

Using these constants in eq. 6, one obtains an ex- 
plicit expression for C(M,M 2 ) in terms of M 

C(M,M 2 ) = C (M, 00 ) = (2-n)rA n-2 M n_1 , r'M £ M 2 (8) 


For impacts between two particles where r'M is greater 
than M 2 , we take the total available mass M + M 2 to be equal M g 

M = M + M 0 (9) 

e 2 


and obtain 


C(M,M 2 ) = (2-n)A n_2 (M + M 2 )M n “ 2 , r'M >_ M 2 (10) 

This relation, together with eq. 2 and 8 defines the model 
crushing law employed in this study. 
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Approximate numerical values for r and A, based on 
hypervelocity impact experiments into basalt by Gault et al 
(1963) are given in Table I at several impact velocities. 


TABLE I 


V( Km/sec ) 


r 



A 


5 

1.3 

X 

10 2 

1.3 

X 

10 

10 

5 

X 

10 2 

5 

X 

10 

1 CT 

1.1 

X 

T A 3 

_L W 

1.1 

X 

1 n 

20 

2 

X 

10 3 

2 

X 

10 


The value of r' is more difficult to estimate. Gault 
(private communication) observed that a basalt particle is com- 
pletely shattered by a projectile 10“ 3 times its mass, moving 
at 2 km/sec. Since r is about 20 at this velocity, 

r' « 50 r 


for this case. 

3 . 0 COLLISIONAL MODEL 


In this section the mathematical formulation of the 
evolution of a system of colliding particles is developed. To 
be specific, given that f(m,t)dm is the number of particles 
having a mass between m and m + dm at a time t, this function 
will change as a result of collisions between the particles 
because many new particles are constantly created and others 
are destroyed. The system itself possesses a "sink" in the 
sense that sufficiently small particles are removed by the 
Poynting Robertson effect and still smaller ones are almost 
instantly blown out of the solar system by radiation pressure. 

In what follows, the system will be assumed suffi- 
ciently random that an effective average collisional velocity 
is meaningful; the collision cross-section is taken as the 
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cross-sectional area of the colliding particles. This assump- 
tion is equivalent to the process of finding the motion of the 
center of mass of the system of particles, then switching to 
the center of mass coordinate system; the particle velocities 
will then be random, to a first approximation. Here we have 
invoked the analogy of a system of gas molecules in a box, 
when the box itself undergoes translation or rotation. 

Assuming spherical particles, the probability of col- 
lision between two particles with radii r^ and r 2 is propor- 
tional to ir(r^ + r^)^. In what follows, the particle masses 

rather than particle radii will be taken as the independent 
variable and whence, the probability of collision per unit time 
between two particles is proportional to 


K 



+ M. 


1/3 


( 11 ) 


where 


K = V 




2/3 


Here V is the average relative velocity of the particles and p 
is the material density of the particles. K will be taken to 
be a constant, to a first approximation. This means that the 
velocity distribution is taken to be independent of the mass 
distribution and all particles are assumed to have the same 
material density. If the expression, eq. 11, is multiplied by 
the number density per unit volume of particles in the mass 
range M 2 to M 2 + dM 2 , the resulting expression 


K 


M. 


1/3 


+ M. 


1 / 3 ' 


f (M 2 ,t)dM 2 


( 11 *) 


is proportional to the total number of collisions (per unit 
time) of an individual particle with mass with other particles 

in the mass range M 2 to M 2 + dM 2 . In other words, the expression 

eq. 11 Is the "influx” per unit time cf particles in the mass 
range M 2 to M 2 + dM 2 "into" a particular object of mass M^. 
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If now the expression, eq. 11 », is multiplied by the 
particle number density per unit volume in the mass range M x 

to lV^ + dM^ the resulting expression 


K 


M. 


1/3 


+ M. 


1/3' 


f(M 2 t)dM 2 


f(M, ,t )dM 


( 11 ’) 


is the total number of collisions, per unit volume and unit 
time, of particles in the mass range M 1 to M 1 + dlV^ with par- 
ticles in the mass range M 2 to M 2 + dM 2 . 


A collision equation defining the collective evolu- 
tion of our system of particles can now be defined. The time 
rate of change of the number of particles in a mass range of 
m to m + dm is given, in a schematic form, by the following 
expression (individual terms are explained below): 


3f (m,t ) 

3 1 


dm 


I 


rate of change of the number of 
particles per unit volume and 
unit time in mass range m to 
m + dm due to erosion 


+ 


II 


+ 


rate of change, because of "ca- 
tastrophic" collisions, of the 
number of particles per unit 
volume and unit time in the mass 
range m to m + dm 


III 


+ 


number of particles in the mass 
range m to m + dm, created per 
unit time and unit volume by 
collisional crushing 


( 12 ) 


IV 


+ 


number of particles in the mass 
range m to m + dm removed per 
unit time and unit volume by 
the Poynting Robertson effect 
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Term I is the rate of change of the number of parti- 
cles per unit volume and unit time in the mass range m to 
m + dm due to the fact that the masses are themselves changing 
in time. This is caused by collisional processes which erode 
particles into and out of the mass range m to m + dm with the 
passage of time. Mathematically, the problem is to find the 
rate of change of the number density of particles in mass range 
m to m + dm at fixed mass , given that the masses themselves 
change at a prescribed rate dm/dt. It can be shown, as is done 
in Appendix A, that the resulting expression for Term I, is 


(I) 


9f (m, t ) 
9 1 


erosion 


9 f (m , t ) dm 
9m dt 


f (m,t ) 


9 dm 
9m dt 


(13) 


For dm/dt, which is the rate at which the mass of a particle 
changes in time, we use the mass removed per unit time by 
collisions with mass M not large enough to completely disin- 
tegrate our mass m. The amount of mass removed during a single 
collision with a mass M is, according to equation 7, 


rM 


(14) 


The number of collisions that one mass m will experi- 
ence (per unit time) with particles in a mass range M to M + dM 
is (cf discussion preceding eq. 11’) 


Kf (M,t ) 




m 


1/3 



(15) 


The total mass removed from m (per unit time) due to 
collisions with particles in the mass range from M to M + dM is 


KTM f (M, t ) 


M 


1/3 



( 16 ) 


Whence, the mass removal rate due to collisions with 
a finite mass range of particles is 
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dm 

dt 



M f(M,t) 



+ m 


1/3 


2 


dM 


(17) 


For the lower limit we take particles with mass y which 
are the smallest particles present. For the upper limit we take 
particles having a mass just large enough to completely break up 
one mass m. As a criterion for "catastrophic" encounter, we take 
(see discussion preceding eq. 8) 


m = r'M 


(18) 


as the limiting mass M still included in the integral (17). 
Therefore, the expression for the mass loss becomes 


dm 

dt 


= - TK 


•m/r 


M 


f(M,t)|M 1/3 + m 1/3 j dM 


(19) 


and term (I) of the collision equation (12) is: 


(I) = 3 " f - 3£ > — TK 


+ FK 


M f(M,t)jM 1/3 + m 1/3 | dM 

C m// ^ 1 I \ 2 

M f(M,t) M 1/3 + m 1/3 dM (20) 


This completes our derivation of term (I) in the collision equa- 
tion. 

We now consider term (II), catastrophic collisions, of 
eq. 12, which can be derived by noting that the probability of 
collision per unit volume and unit time between two particles in 
a mass range m to m + dm and M to M + dM is (cf discussion 
accompanying eq. 11’') 


m 


1/3 


+ M 


1/3 


f(m,t) dm f (M,t )dM K 


2 


( 21 ) 
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The probability for a collision per unit volume and 
unit time between a particle with mass m and some other parti- 
cle In a finite mass range is then 


2 

Kf (m,t ) j f(M,t)|m 1/3 + M 1/3 J 


dM 


( 22 ) 


Since the effect of collisions with particles that do 
not completely destroy m has been accounted for in term (I), we 
evaluate the integral in eq. 22 over all "catastrophic” colli- 
sions. The result is term (II): 


r 


(II) = - Kf (m, t ) 


M 


f(M,t) 


m/r 


m 1/3 + M 1/3 dM (23) 


where the minus sign is used to denote a particle removal pro- 
cess. 


Piotrowsky (1953), in an earlier study, obtained a 
term similar to eq. 23. We disagree with him, however, inas- 
much as Piotrowsky' s formulation is equivalent to 

r M w 

-Kf (m, t ) m 2/3 f (M,t ) dM 

Jm/I” 

1/3 1/3 2 

i.e., the collision cross-section factor (m + M ) is re- 
placed by m 2/ ^ 3 . This approximation is invalid for M > m, and 
for arbitrary distributions can introduce a serious error near 
the upper limit. 

Term (III) of the collision equation (eq. 12) can be 
derived by noting that for each event whose probability is given 
by eq. 22 there are g(m; M,M 2 )dm "secondary" particles produced 

in the mass range between m and m + dm, where g(m;M,M 2 ) is the 

crushing law eq. 3 and where M is the smaller one of the two 
colliding masses and is assumed to be completely broken up 
during impact. One, therefore, has 
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dmK 


= Km" n dm / dM / dM 2 C(M,M 2 ) f(M,t) f(M 2 t) 


I dM J~ dM 2 g(m;M,M 2 ) f(M,t) f(M 2 t)|M 1/3 + M 2 1/3 


M 1/3 + M 2 1/3 


(24) 

\ 2 


which is the rate at which particles having a mass between m and 
m + dm are created (per unit volume and unit time) due to the 
crushing of mass during an impact between sufficiently large 
particles to produce fragments of mass m, and where C(M,M 2 ) is 


given by 


firtn t nn o K r»r» 
x^iiu G L/l 


i n 

J- V I 


The limits on the double integral in equation 24 can 
be obtained by noting that the dummy M is always smaller than 
dummy M 2 , otherwise, the argument of C(M,M 2 ) should be replaced 

by M 2 ,M (this follows from the definition of the crushing law, 

discussed in the previous section) . The lower limit of the 
integral for M 2 is, therefore, M. Since we are interested in 

all possible collisions, the upper limit for M 2 is M o . The 

lower limit for M is defined by the condition that the largest 
fragment M^ produced during impact must be equal to or larger 

than m whence, by equation 7. 


m. < M, = AM 
— b 


(25) 


and 


M > m/A 


( 26 ) 


Using appropriate upper limits and dropping the dif- 
ferential dm, equation 24 becomes term (III) in the defining 
equation 12: 
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(III) = K(2-n)m n A n ~ 2 A 2 x 


-M^/r ' ^r’M 

x< | dM dM, 

/m/A Jm 


(M+M 2 )|m 1/3 + M 2 1/3 | M T1 “ 2 f(M,t)f (M 2 ,t ) 


-M /I” r M_ 


+ r dM 

m/A 


dM, 


r»M 


M 1/3 + M 2 1/3 J M n " 1 f(M,t)f (M 2 ,t) 


,-M 


1 

00 

dM 

oo 

dM 2 (M+M 2 ) 

. 2 | 

+ M, 1/3 

' 

^ M /r ’ 
00 

M 



(27) 


where the first and third integrals refer to catastrophic colli- 
sions between masses M and M p such that both are totally disrup- 
ted 


r 1 m >_ m 2 


( 28 ) 


the second integral refers to erosive collisions between masses 
M and M 2 such that M p behaves as an infinite target 


r'M £ m 2 


and the mass redistributed is just rM. 

The third integral refers to catastrophic collisions 
between objects in the mass range M^/r* to M^. This means that 

for test masses m 

m/A > M^/r ' 


the first two integrals are zero and only the third integral is 
retained* with lower limit of m/A replacing M^/I”. 


*See Appendix B, for details. 
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The last term (i.e., term (IV)) In eq. 12 Is the rate 
at which particles are being lost due to the Poyntlng Robertson 
effect (Robertson, 1937) due to radiation damping. According 
to Robertson's analysis, a particle in an approximately circular 
orbit will move from a mean distance from the sun of R q (AU) to 

a distance of R(AU) from the sun, during a time t given by 


t 




R 2 ) m 1/3 


years 


(29) 


where a particle specific gravity of 3*5 has been assumed, and 
where m is expressed in Kilograms. 

For an effective particle lifetime one may take the 
time required for moving 1AU closer to the sun (R = R - 1 

although this is arbitrary to some extent), in which case, one 
has 


t = 10® (2R q - 1) years 

If R = 1AU, then x becomes 
o 

t = 10® m 1 ^ years = t q m 1 ^ (30) 


The rate at which the population per unit volume of 
particles having a mass between m and m + dm decreases in time 
can then be represented approximately as 


(IV) 


f (m,t) 

T 


-m '*~ / '^f(m,t) 


(3D 


where use has been made of eq. 30, and where the differential 
dm has been dropped. 

4.0 SOLUTION OF THE COLLISION EQUATION FOR SMALL PARTICLES 


The collision equation derived in the previous sec- 
tion is complicated. Simplifications can, however, be obtained 
by restricting one's attention to special cases. We shall, in 
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what follows, consider special cases where the collision equation 
can be solved. We will limit our attention to the evolution of 
particles having a mass m very much smaller than M^, but very 

much larger than y which is the smallest particle mass present 
and is determined by the radiation pressure limit. 

More specifically, particles with a mass smaller than 
r’y will no longer be eroded according to this model, since any 
collision they experience will be ’catastrophic’. Our present 
model is, therefore, expected to be valid only for masses m in 
the range 


r’y << m << AM /r ' 

o® 


(32) 


where r’y is the largest mass completely disrupted by a colli- 
sion with the smallest surviving particle y; this means that a 
mass <_ r’y is completely disrupted by every collision and the 
erosion term for particles in this mass range does not apply. 
AM^/r' is the largest mass that can be created by a catastrophic 

collision and therefore the particle creation term (III) does 
not apply to the distribution of masses greater than AM^/r'. 

The collision eq. 12 expresses the time rate of change 
of the particle number density function in terms of the individ- 
ual collision processes. We shall choose the simplest possible 
case and as a first approximation assume that the particle dis- 
tribution has reached "steady-state", i.e., the rate at which 
particles having a mass m are produced is balanced by the rate 
of their removal by collisions. 

Mathematically, the steady-state assumption implies 
that a source exists which steadily regenerates the large objects. 
This appears to be the case with cometary meteoroids, since comets 
are believed to be constantly giving off particles. While the 
rate of generating particles by comets as a function of particle 
size is not precisely known, one has some information on the sub- 
ject by considering stream meteors which are believed to have been 
generated comparatively recently. 

Radio work indicates that the abundance of faint 
meteors in streams is small by comparison with the background 
sporadic flux, while nearly half of the bright photographic 
meteors move in identifiable streams (Whipple, 1963). It, 
therefore, appears that the meteoroids generated by comets are 
comparatively "rich" in large particles and poor in small ones. 
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We shall, to a first approximation, assume that the only effect 
of the meteoroid production by comets on our model distribution 
is one of keeping fCM^) constant in time, i.e., keeping the 

"top masses" of the distribution constant. 


In the case of asteroids, the situation is more specu- 
lative since large asteroids are rot believed to be regenerated 
at the present time. If, however, we assume that the processes 
of particle creation by fragmentation and particle removal by 
collisions are approximately equal then the present theory is a 
reasonable first approximation. Among the catalogued asteroids, 
certain groups exist (Hirayama groups) that are believed to have 
resulted from collisions by "parent" asteroids. Such collisions 
between very large asteroids are believed rare, with a lifetime 
of perhaps many millions of years. As the size of the asteroids 
diminishes, however, their frequency increases with a correspond- 
ing increase in the number of collisions. 


We now turn our attention to the individual terms of 
the collision equation (eq. 12). Their properties will be dis- 
cussed in detail and appropriate simplifications will be intro- 
duced. 


Term (I) (eq. 20) of the collision equation expresses 
the rate of change due to erosion in the particle number density 
function and has the mathematical form: 


(I) = 


8f (m,t ) K 

3m 


r m/T' 


Mf (M,t ) 


M 1/3 + m 1/3 | dM 


+ f(m,t)TK 


m 


( r ') 2 


Sr.tl <rT 1/3 


+ 1 


m 2/3 (33) 


+ f (m, t ) rK | m -2/3 


m/r 


Mf (M,t) 


M 1/3 + m 1/3 j 


dM 


where the last two terms are the result of carrying out the par- 
tial differentiation with respect to m, as indicated in eq. 20. 
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Since r* is of the order 
sional cross-sectional area can be 


1 4 

of 10 J to 10 , 

taken to be 


the colli- 


K 


M 


1/3 + m 1/3 


^ Km 


2/3 


(34) 


The error thereby committed in the upper limit is, 
(^,1/3 + m l/3) 2 . m 2/3 . m 2/3 (2r -l/3) % 10% of m 2/3 (35) 

and the error committed in the lower limit is completely negii- 
ble. We therefore write, approximately 


(I) £ i f 0q>*) .rK m 2/3 

din 


m/r’ 


Mf (M,t)dM 


+ f(m,t) f(m/r’,t) rK/(r ’) 2 m 5/3 


(36) 


+ f(m,t)TK (2/3) m 1/3 


rffl/r' 


Mf (M,t)dM 


where the above approximation is a slight "under estimation" of 
the process to the extent that the geometrical cross-section of 
the smaller particle is neglected. The approximation is, how- 
ever, far less serious than may appear on the surface, because 
the cross-section 


K 



+ m 


1/3 


2 


includes the case of "grazing incidence"; and, therefore, over- 
estimates the probability of mass removal during cratering. We 
are, therefore, satisfied that equation (36) defines process (I) 
with reasonable accuracy. 
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We now consider the rate of "catastrophic” collisions 
given by term II of the collision equation (eq. 23): 


f M oo 

-Kf (m,t ) I f(M,t) 

J m/r 1 


m 1/3 + M 1/3 dM (23) 


Using a power law steady-state distribution 


f(m,t) = f(m) = Am “ 


( 1 ) 


and eq. 23 integrates to give (multiplying out the squared 
parenthesis) : 


- KA 2 m” a jm 2/3 


M " a+1 - (r ’) 01 " 1 m -a+1 


( -a+1) 


-1 


+ 2m 1 / 3 [m -a + V3. (rl) a-V3 m -a + V3 


(-a+4/3) 1 (23*) 


+ [ M co " a+5/3 - (r*) a " 5 / V a+5/3 


(-a+5/3) 




AM 


provided that a 4 5/3, 4/3 or 1. If now a > 5/3 and ^r 
this expression simplifies to 


m 


- KA 2 m -a m “ a+5/3 


x constant 


and if 


a < 5/3 


AM 

00 

r'm 


oo 


one obtains, instead. 


- KA 2 m-“ M -“ +5/3 


x constant. 
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It therefore appears that for a population of a power 
law type the collisional lifetime of an object is only a func- 
tion of the more numerous small objects for a population index 
a > 5/3 and for a population index of a < 5/3 it is given by 
the -a + 5/3 power of the largest object in the sample. 

The collision equation (eq. 12) can be further simpli- 
fied when the relative importance of radiation damping (term IV) 
is considered. We shall show, in what follows, that this term 
can be neglected altogether by comparison with the influence of 
collisional processes on the particle number density function. 
The lifetime of a particle limited by the Poynting Robertson 
effect is, according to equation 30, 


in 8 1/3 

x pR =10 m year 


(38) 


where the subscript PR stands for Poynting Robertson and where 
m is the particle mass in Kilograms. 

The lifetime of a particle of mass m due to catastro- 
phic collisions by another particle with mass m/r' or larger, 
can be estimated by noting that 


T CC 


1 

Na 


(39) 


where the subscript CC stands for catastrophic collisions, N is 
the flux of particles having a mass of m/r’ Kg or larger per 

meter 2 sec and a is the cross-sectional area of the particle in 

meter 2 . Taking r* £ 10** , particle density of 10 3 kg/meter 3 and 
the NASA model* cometary flux of 


N 



6 



(40) 


one obtains 


T CC 


£ 2 x 10 6 m 1/3 


years 


where m is the mass in Kg. This is 100 times shorter than 

t PR' 


* Natural Environment and Physical Standards for the Apollo 
Program, April 1965 - revised. 
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Another check on the relative importance of radiation 
damping may be performed by estimating the lifetime of a parti- 
cle T p limited by erosion. The mass erosion rate of a cometary 

particle has been estimated (Whipple, 1963) as 7 x 10 

2 

gm/cm sec. This means that we may take, in MKS units. 


dm 

dt 


= - b x 


4irr = - 


b x 5 x 10 2 m 2/3 


(41) 


where dm/dt is the rate of mass loss due to etching, b is the 
rate of mass loss per unit surface area and where a particle 

density of 10 J Kg/m J (1 gm/cm- 3 ) has been assumed. 

This equation is easily solved for the time re- 
quired to erode the particle: 


t E ^ 


3 x m 


1/3 


5b x 10 


-2 


(42) 


Using b = 7 x 10 


-12 


Kg/meter sec 


we get 


t E - 


10 


13 


sec = 3 x 


10^ m 1//3 year 


(43) 


again yielding a process two orders of magnitude faster than 
radiation damping. These results are in agreement with an 
earlier discussion by Whipple ( 1963 ) where he treated the life- 
time of photographic meteors. 

We are, therefore, justified in disregarding term IV 
in the collision equation for particle distributions near earth. 
In order to avoid the introduction of special and arbitrary 
assumptions, the status quo of the meteoroid environment in the 
near earth space will also be extended into the asteroidal belt 
under the present model; this means that the lifetime of the 
asteroidal debris is taken to be their collisional lifetime 
rather than their lifetime due to radiation damping. Term IV 
in the collision equation is thereby disregarded altogether. 
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Using the steady-state assumption and equations 36 , 

23 and 27, substituted into equation 12, the collision equation 
can now be expressed as: 


af(rc.t) , „ _ v „2/3 „ 9f(m) 


St = 0 = r -jg- 


.m/r' 


Mf (M)dM 


+ (2/3) f(m)rK m ~ 1/3 


m/r’ 


Mf (M)dM 


+ f(m) r(r') 2 K f(m/r') m 5/3 

r M » 

- K f (m) I f(M)jm 1/3 + M 1/3 ] 

Jm/r* 

+ K(2-n) m~ n A n ~ 2 x< 

-M /r* 


dM 


(4K) 


M oo /r ’ r r ’ M 

dM / dM 2 Q(M,M 2 ,m) 
m/A J M 


+ r 


dM 


m/A 

rM 


dM 2 R(M,M 2 ,m) 


r’M 


M 


dM dM 2 Q(M,M 2 ,m) 

M /r* J 


where the time dependence has been dropped from f(m,t) and where 

Q(M,M 2 ,m) = (M+M 2 )Jm 1/3 + M 2 1/3 j M n “ 2 f(M)f(M 2 ) 

(45) 

R(M,M 2 ,m) = |m 1/3 + M 2 1/3 j M n-1 f(M)f(M 2 ) 
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Seeking a simple power law solution of the form m 

we take 


f(m) = A m“ a 


(46) 


AM 

CD 



and substitute it into equation 12. 

Performing the integrations, one obtains. 


0 = KrA^ 1 (2/3-ot) -a-1/3 

2 - a 


(m/r') 2 ' a - ii 2 ‘° 


+ Kr(r ')“" 2 A 2 nf 2a+5/3 


(47) 


KA 2 (r ') a ~ 1 m -2a+5/3 

a - 1 


+ KA 2 C* m " 2a+5/3 


provided that both of the following inequalities are satisfied: 


a > 5/3 


a > ± (n + 5/3) 


(48) 

(49) 


The quantity C' is a lengthy expression involving a, n, r, r* 
and A. 


If a < 5/3 > no solution of a population index type 
exists, as can easily be shown; the dynamic processes of parti- 
cle creation and catastrophic collision then depend on the 
largest masses present in the distribution and the linearity 

of eq. 47 in m -2a+ 5/3 is destroyed. A similar thing happens 

when a > 2; the term y in the first term of eq. 47 will 
dominate because y << m/r' and we obtain, for this erosion pro- 
cess 


const x m “ (0 - y 2 “) 


(50) 
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Equating each power of m in eq. 47, one obtains 


-2a+5/3 = - a - 1/3 


resulting in 


a = 2 


yielding the absurd result that if a > 2, then a = 2, indicating 
that for a > 2 no solution of a population index type exists. 

Experimental work by Gault et al ( 1963 ) on hyper- 
velocity cratering into basalt indicates that n % 1.8. Some 
cases with slightly higher n (but still less than 2) have been 
observed together with lower values for n . For very small 
particles, the value for n appears to be still lower (Gault et 
al 1965 * Gault and Heitowit, 1963 ). While these experiments 
refer to basalt only, the writer is not aware of any other 
determination of the size distribution for crushed debris 
during hypervelocity cratering. It is a significant but not 
necessarily serious mathematical limitation of the present 
model that if the parameter n is equal to or is larger than two, 
the form of the crushed particle distribution changes, as can 
be seen from equation 12 and 13 when n > 2 is substituted. In 
that case, the quantity C(M,“) becomes, approximately 


C(M,») = 


(n-2)M e 

" 7 ^ 


(n-2) rM 
2-n 


(51) 


and the discussion in this paper is no longer applicable. 

We now consider, briefly, the time dependent colli- 
sion equation (eq. 44, with 3f/3t 4 0)* When a separable popu- 
lation index type of solution 


f(m,t) = A(t) m “ 


is substituted into eq. 44, we obtain (with the help of eq. 47 ) 
m~ a = constant x A 2 (t) x m - 2 a+ 5/3 
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This equation is satisfied only if a = 5/3. It can, however, 
be shown readily that the particle creation term III (eq. 27), 
in this case becomes (after iteration with a = 5/3 and n > 5/3). 


p __ [T / O _r /n 

A (t) x (constant x m -n + terms in m 3 ^ and m ~ >/ D x £n m) 


and the collision equation is obviously not satisfied. We 
therefore conclude that no separable population index type of 
function exists (other than zero) which satisfies the time 
dependent collision equation employed in the present study. 

It is therefore concluded that within the present 
model, a population index type solution exists if and only if 
the population has reached a steady-state distribution and that 


r ' u << m 


AM 

00 



(52) 


\ (n+5/3) % 1.75 < a < 2 


in which case eq. 47 reduces to 


0 


constant x m 


-2a+5/3 


(53) 


In this equation, the constant factor is an algebraic as well 
as transcendental expression in a and the physical parameters. 
Knowledge of the latter then permits one to calculate the 
former. Such a study together with a discussion of the sta- 
bility of the solution is presently under preparation. 

An interesting property of the solution in the mass 
range specified by eq. 52 is that, the particular mass ranges 
which dominate dynamically the number density of objects in 
the mass range m to m + dm are: 
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Term (I): (Erosion) is "dominated" by objects of mass m/r* 

Term (II): (Catastrophic collision) is "dominated" by ob- 
jects of mass m/r' (53) 

Term (III): (Particle creation) is "dominated" by masses of 
the size m/A impacting masses of the size r’m. 


In other words, the dynamics of the collisional pro- 
cess described by the present model is to a first approximation 
independent from the size of the "cutoff" objects y and M^. 


This is an important conclusion. It indicates that slow deple- 
tion of the largest objects (e.g., asteroids) will not invali- 
date the collisional model in a first order (of approximation) 

r'PMt) 


for masses within the range r'n << m << 


where M (t) 


designates the mass of the largest objects present as a function 
of time. 


5.0 APPLICATION TO THE NEAR EARTH METEOROID ENVIRONMENT 


The results of this paper are compared with experimen- 
tal information in Figure 1. The figure is a plot of the ex- 
ponent a as obtained from the present model 


1.7 < a < 2 


( 5*0 


versus particle mass. 

The small mass limit of the model is 


m > > r ’ y 


(55) 


Gault (private communication) observed that a basalt particle 

is broken up catastrophically by a projectile 10 times its 
mass moving at 2km/sec. Since meteoroids move with a relative 
velocity of about 15 to 20 km/sec and since cometary particles 
are believed to break up more easily than basalt, we take ten- 
tatively 


r ' 


<\/ 


10 


5 


(56) 
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■3 3 

The limiting mass with a density of 10 kg/m not 

-15 

blown away by radiation pressure is about 10 kg. We there- 
fore have. 


m >> 10 10 kg (57) 

This means that for masses 

m < 10~ 10 kg ( 58) 


the present model is not applicable. 

It is interesting to note that in precisely this range 
of masses a definite change in the distribution is observed 
(Figure 1). The "slope" of the logarithmic distribution changes 
to smaller absolute values. Near the radiation pressure limit, 
the slope should move toward a positive value, indicating that 
the number density of particles increases with increasing mass 
and at the radiation pressure limit we expect a singularity 
(<*-»•“>) reflecting the fact that no particle should exist having 
a mass equal to the radiation pressure limit or smaller. 

For masses 


m > 10 


-9 


kg 


the radar results are a lower bound in the range of values for 
-a obtained in this model. The value of a = 2.3^ suggested by 
Hawkins and Upton (1958) is somewhat high, but another analysis 
(Dohnanyi, 1966 and 1967 ) of comparable photographic meteor 
data (McCrosky and Posen, 1961) indicates a slope a = 2, in 
agreement with the present model and irons (Hawkins, i 960 ) have 
a = 1.5 which is somewhat low. 

It is significant that the range of values for a in 
the present model includes the experimental value of the popu- 
lation index as its upper limit. It is strongly implied that 
cometary meteoroids undergo the collisional processes discussed 
in this paper and our model is a first approximation to their 
distribution. 
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6.0 APPLICATION TO ASTEROIDAL DEBRIS 


This section is a comparison of the present results 
with the distribution of the known asteroids. While a detailed 
comparison requires a discussion of the time dependent problem 
(because the largest asteroids are no longer being created) 
which is beyond the scope of the present study, it is interest- 
ing to note the comparison between the present results and 
observation. 

Figure 2 is a logarithmic plot of the cumulative 
number of catalogued asteroids as given by Kuiper et al (1958) 
versus absolute photographic magnitude. An albedo of .1 and a 

1 3 

density of 3-5 x 10 Kgm/m J have been assumed in associating an 
asteroidal mass with a given absolute photographic magnitude g= 
The following relationship can then easily be derived 


m 


10 


24- . 6g 


(59) 


where m is in Kg. 

The straight line, in this figure, is the result of 
a least squares fit to the unnormalized distribution where the 
3 largest asteroids have been disregarded. The population index 
a of the fit is 


a = 1 . 80 + .04 


( 60 ) 


The dashed line indicates the "true number" of asteroids esti- 
mated* by Kuiper et al on the basis of their analysis of selec- 
tion effects influencing the probability of detection. The 
reason why the distribution approaches a horizontal line for 
sufficiently faint asteroids is believed to be due to such an 
effect . 


*The "true number" of asteroids has been estimated by 
Kuiper et al over a much larger range of asteroids (to an 
absolute magnitude of g = 13), but in view of uncertainties, 
we only consider the estimated distribution for g <_ 11. 
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It can be seen, from the figure, that a straight line 
is a good fit to the data in the region 


6 < g < H (61) 


if we correct the distribution for selection effects in the 
region of 10 _< g _< 11 as suggested by Kuiper et al. The fit 
is good, even if data in the region 


6 <_ g <_ 9*5 


( 62 ) 


are only used. As seen from Fig. 2 the three largest asteroids 
deviate from the trend exhibited by the 647 others considered. 
Since for large asteroids no sources of generation are known to 
exist, the steady-state assumption could hold only for masses 
considerably smaller than the largest ones present. Thus, the 
distribution of the known asteroids appears to be consistent 
with the results of the present model. 

The distribution of the lunar craters* yields, on the 
average, a population index 


1.6 < a < 1.7 


(63) 


which is close to but a trifle lower than the result 


1.75 < ct < 2 


(64) 


obtained in this study. 

Hawkins (I960) has studied the distribution of meteo- 
rites and obtained 


a = 2 and a = 1 . 5 


( 65 ) 


*See, for example. Fielder (1963), Dodd et al ( 1963 ), 
Brinkman ( 1966 ) Hartmann (1964) and Baldwin (1964). 


Bellcomm, Inc. 


- 29 - 


depending on the kind of meteorite (i.e., stones or irons, re- 
spectively). Stones are therefore seen to have a population 
index in agreement with the present model and irons have a some- 
what lower population index. 

7.0 CONCLUSION 


A collisional model of interplanetary debris is formu- 
lated and solutions of a population index type are sought, i.e., 
a solution for the number density of the time independent form 


f(m)dm = m a dm (66) 

is substituted into the collision equation (eq. 12). 

It is found that eq. 66 indeed solves the collision 
equation in the mass range of 

r’y << m << AM^/T ’ ( 67 ) 


provided that the population has reached a steady-state dis- 
tribution. The quantity y is the largest object completely 
shattered by a particle of mass r and AM^/r' is the largest 

fragment produced when an object having a mass is "cata- 
strophically" shattered. 

The population index a must satisfy 


| (n+5/3) = 1-75 < a < 2 ; (68) 

beyond the range of a, the collision equation (eq. 44) has no 
population index type solutions. For a > 2 erosion by the 
smallest particles dominates and depletes the mass distribution 

with time. For a < ^ (n+5/3) collision products from the most 

massive objects cause evolution of the mass distribution with 
time. It is furthermore shown that, if the population is 
evolving rapidly so that 3f(m,t)/9t 4 0 in e< 4* 44, no separable 
population index type of solution exists, that would satisfy 
the time dependent collision equation. 
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Application of the present results to cometary mete- 
oroids yields satisfactory agreement in both a and the trend of 
the distribution near the small mass end of cometary meteoroids. 
This is indicated in Figure 1. We note that the population in- 
dex for these small particles obtained experimentally is in the 
neighborhood of 


a = 2 


(69) 


i.e., at the upper limit of a as obtained in this study. 

A least squares fit to the distribution of asteroids 
catalogued by Kuiper et al (1958) yields (Fig. 2) 


a = 1.8 + .04 


(70) 


which is within the range of values for a permitted by this 
model. 


Lunar craters, having a distribution with a popula- 
tion index in the approximate range 


1.6 < a < 1.7 


are a trifle below the lower limit of a obtained in this study. 


The result of this model are in reasonable agreement 
with the distribution of stony meteorites as obtained by 
Hawkins (a = 2) and the iron meteorites (a = 1.5) have a rather 
low population index. 



J . S . Dohnanyi 
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GREEK CHARACTERS 


a = population index. 

r 1 = threshold mass of finite target completely 
shattered, per unit projectile mass. 

m 2 = T'M 


where 


M 

“2 


is the target mass and M is the projectile 
mass . 

r = total ejected mass from a semi-infinite target 
per unit mass of the incident particle. 


M e = TM (eq. 7) 


where M is the mass of the incident particle. 

A = mass of largest fragment per unit mass of 
incident particle. 


M b = AM (eq. 7) 

n = population index of ejected particles, (see 
eq. 2) . 

y = mass of smallest object in the solar system 
not blown away by radiation pressure. 

p = material density of the particles. 

t = particle lifetime. 
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LIST OF SYMBOLS 


a 

b 

c(m,m 2 ) 

f (m,t )dm 


g 

g(m;M,M 2 )dm 

k 


m,M 



r 

R 

t 

V 


2 

particle cross-sectional area (m ) . 

coefficient, defined by eq. 4l. 

normalization coefficient of the comminu- 
tion law. 

particle number density function, i.e., 
number of particles in the interval m to 
m + dm per unit volume of real space at 
a t ime t . 

absolute photographic magnitude. 

comminution law for ejecta (eq. 2). 

parameter, equal to • 

particle mass (in Kg). 

mass of the largest ejected fragment. 

total ejected mass during cratering by 
impact . 

mass of largest object in the distribution. 

particle radius. 

distance from the sun in AU. 

time . 

average relative velocity of colliding 
particles . 
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APPENDIX A 


The purpose here is to derive an expression for the 
partial derivative with respect to time of the particle number 
density function due to the process of erosion (i.e., term I 
in the collision equation). More specifically, given that: 


f (m,t )dm 


( A-l) 


is the number of particles having a mass in the range m to 
m + dm at a time t, and that the masses change with time 
according to the relation 


dm 

dt 


€(m,t ) 


( A-2 ) 


The problem is to find 


lim f(m,t+dt) - f(m,t) _ 9f (m,t ) 
dt+0 dt 9t 


due to the process defined by eq. A-2. 

In order to obtain an expression for eq. A-3 we con- 
sider the following equalities: 


f (m,t ) Am = number density of particles per unit (A-l]) 

volume in the mass range m to m + Am ' 

f (m,t ) 5 (m,t )dt = number of particles (per unit 

volume) which have left the mass range (A-5) 
m to m + Am during a time dt 

-f (m+Am,t)€(m+Am,t)dt = number of particles (per unit 

volume) which have "entered" the 
mass range m to m + Am during dt 
provided that Am>£(m+Am,t )dt=dm. (A-6) 

Since, however. Am is arbitrary, 
we choose Am such that it satis- 
fies this condition. 
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Appendix A (cont'd) 

Using eq. A-3, 5 and 6, we have 


9f(m,t) A „ _ lim -f (ra+Am,t ) 5 (m+Am,t )dt+f (m,t ) £ (m,t )dt , „ 
9t ~ dt+0 dt 


Since Am can be taken to be small, we can write 


f(m+Am,t) = f(m,t) + ■ — Am 

o m 


( A-8 ) 


and 


5 (m+Am) = £(m,t) + — Am 


(A-9) 


Substituting eq's A-8 and 9 into A-7 one has, after 

— 2 

expanding the product and dropping the term in Am 


9 f (m, t ) _ -9_ 
9t 9m 


f (m,t)£(m,t ) 


-9 

9m 


f (m,t ) 


dm 

dt 


( A-10) 


which is the desired result. 
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APPENDIX B 


This appendix is a discussion of some of the proper- 
ties of the integrals occurring in term III (eq. 27) of the 
collision equation. 

First we derive the limits of integration of the 
expression eq. 24: 


where 



P(M,M 2 ,m) 


(B-l) 


P(M,M 2 ,m) = Q(M ,Mg ,m) 

when M > Mg/r ’ 

= R(M,M 2 ,m) 

when M < Mg/r ’ 


with Q and R given by eq. 46 in the text. 

We first note that we want to integrate over all M, 

and M 2 for 


M <_ Mg 


(B-3) 


such that particles of mass m are produced. 
This means that 


M < M 0 < M (B-4) 

— 2 — 00 

and 

? < M < M (B-5) 

A — — oo 
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Appendix B (cont'd) 
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Since, however, the value of the kernel under the 
integral eq. B-l has the dependence indicated in eq. B-2, some 
additional restrictions exist. The situation is indicated in 
Fig. B-l, which is a plot of the regions of integration in M 
and M 2 • Dashed lines through the origin are the loci of points 

where 


m 2 = r'M 


and 


M 


“2 


s M 


(b-6) 


(B-7) 


These lines separate the M 2 x M space into regions (with the 
bounds on M and M 2 indicated) with the following significance: 


Region 

1, 

M 2 

< r'M 

II 

<0 


Region 

2, 

M 2 

> r'M 

P = R 

(B-8) 

Region 

3, 

m 2 

< r'M 

►n 

n 

<0 



whence, integrating over the various regions, the integral eq. 
B-l becomes 


• M /r' /-r’M 


K m -n / dM I dM 2 P = K m -T1 


• M /I" 


+ r 


dM 


Wa Jr'M 


■M 


dM 2 R + 


‘m/A 
•M 

CO 

Jm /r 


dM 


dM 2 Q 


dM 


'M 


•M 


(B-9) 


dM 2 Q) 


M 


provided that 


m/A < M /r» 


(B-10) 


as can be seen from the figure. 
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Eq. B-9 overestimates the mass production process at 
the lower limit in region 1 inasmuch we have included processes 
where two objects, both having a mass of m/A, will create frag- 
ments of the size m which is an absurd result. A more detailed 
treatment would consider the precise expression for the mass of 
the largest fragment when two objects of similar mass collide. 





I 


FIGURE I 












FIGURE 2 
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